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Abstract 

We investigate the statistical properties of the scattering matrix S describing 
the electron transport through quasi-one dimensional disordered systems. For 
weak disorder (metallic regime), the energy dependence of the phase shifts of 
S is found to yield the same universal parametric correlations as those char- 
acterizing chaotic Hamiltonian eigenvalues driven by an external parameter. 
This is analyzed within a Brownian motion model for S, which is directly 
related to the distribution of the Wigner-Smith time delay matrix. For large 
disorder (localized regime), transport is dominated by resonant tunneling and 
the universal behavior disappears. A model based on a simplified description 
of the localized wave functions qualitatively explains our numerical results. In 
the insulator, the parametric correlation of the phase shift velocities follows 
the energy-dependent autocorrelator of the Wigner time. The Wigner time 
and the conductance are correlated in the metal and in the insulator. 

PACS numbers: 72.15.-v, 73.20.Dx, 72.10.Bg, 05.60,+w 
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I. INTRODUCTION 



A. Preface 



Electron transport through quasi-one dimensional disordered systems is characterized 
by the scattering matrix S, relating the amplitudes of incoming and outgoing waves. The 
scattering matrix contains information about not only quantities as the conductance and 
the characteristic dwelling times for electrons moving through the disordered region, but 
also on the energy levels of the system. The universal statistical properties found in the 
energy spectrum of disordered and chaotic systems have been related to the distribution 
of eigenvalues of random matrix Hamiltonians On the other hand, the universal 

conductance fluctuations of metallic systems have been understood in terms of a statistical 
description of the transfer matrix [§,0. To understand the relation between these two 
universalities, it is useful to study the statistical properties of the scattering matrix HBHTT[|. 

Due to current conservation, the scattering matrix is unitary and its eigenvalues are rep- 
resented by 2N phase shifts {0{\. (N is the number of transverse channels in each asymptotic 
region.) Assuming that all matrices S with a given symmetry are equally probable, we ob- 
tain Dyson's circular ensembles fT2] , |T3|| (named COE and CUE for orthogonal and unitary 
symmetry classes, respectively), where the phase shift statistics follows universal laws. The 
isotropy hypothesis of Dyson's ensembles applies only to ballistic chaotic cavities with no 
direct channels |14,8||, where the electronic motion is essentially zero- dimensional after a 
(short) time of flight. The conductance in this systems is always of the order of N/2 since 
reflection and transmission are equally likely. On the other hand, in disordered systems this 
isotropy hypothesis is not satisfied any longer because transmission is much less probable 
than reflection. For quasi-one- dimensional (quasi-ld) metallic samples the conductance is 
of the order of Nl/L, where / is the elastic mean-free-path, L the length of the disordered 
region, and I L. In the localized regime, when the localization length £ ~ Nl is much 
smaller than the sample length L, the typical conductance scales as exp (— 2L/£). However, 
the failure of Dyson's hypothesis for the mean values of S (which is linked to the average 
conductance) does not prevent the fluctuations of S for weakly-disordered quasi-ld samples 
to be (approximately) described with the universal behavior of the circular ensembles, as 
shown in Ref. ||11|| . This good agreement gets poorer if the disorder is increased, as well 
as when the system does not have a quasi-ld geometry, or when it enters the localized 
regime. In the first two cases the mean values of the transmission and reflection amplitudes 
become strongly dependent on the channel index and the eigenphase distribution is highly 
anisotropic. When localization is achieved by keeping the disorder weak and increasing the 
length of the sample beyond the localization length, the scattering matrix was shown to 



decouple in two statistically independent (almost unitary) reflection matrices jy]]. This 



picture is particularly clear if one uses a semiclassical approach where most of the classical 
trajectories return to the region of departure instead of traversing the disordered sample. 

A few years ago, another type of universality in the spectrum of chaotic and disordered 
systems was discovered by Szafer, Altshuler, and Simons |T3|,[TB|. It concerns the adiabatic 
response to an external perturbation (magnetic field, shape of the confining potential, etc). 
The correlator of the derivatives of the eigenenergies at two different values of the external 
parameter has a universal functional form once a proper rescaling is carried out. Latter 
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studies have found these universal parametric correlations to hold in a variety of other 



systems, including interacting one- dimensional models ||17|| , ^-dependent band structures 
of semiconductors |18[], Rydberg atoms in magnetic fields |]T9fl , etc. Considering the 2N 
eigenphases of S as a function of the electron energy, the natural question arises as to whether 
the analogous parametric correlations will have the universal form found for eigenenergies. 
In the weakly-disordered quasi-ld case, where the fluctuations of the eigenphases follow 
the circular ensemble predictions and N ^> 1, one could anticipate that the parametric 
correlations do exhibit the universal form of Szafer et ai, and this is confirmed here by 
numerical simulations and analytical arguments. More interesting is the question of what 
happens to these correlations as we go into the localized regime, and whether or not they 
can be described by a simple model as that of the decoupling of the eigenphases. 

Directly related to the parametric correlations of the eigenphases, there is the question of 
the statistical distribution of the traversing-times in the disordered region. This distribution 
and the related correlation functions have recently begun to be addressed in Refs. PU||?T|| , 
because the Wigner time appears in a variety of physical phenomena, ranging from the ca- 



pacitance of mesoscopic quantum dots |22j to nuclear resonances [^3j . Since the Wigner time 
is obtained from the scattering matrix and its energy derivative, the statistical properties 
of S(E) determine the correlations of the traversal time. We aim in this paper to study the 
relationship between the parametric correlations and the distributions of the traversal times, 
and their connection to the statistical properties of the corresponding scattering matrix for 
metallic and localized quasi-ld disordered systems. 

In the remaining of this section we introduce the basic notation concerning the scattering 
matrix and present the substantially different behavior of the eigenphases, conductance, and 
Wigner time (as functions of the Fermi energy) between the metallic and localized cases. 
In Sec. [TJ we study the parametric correlations of weakly-disordered quasi-ld systems. We 
justify their universal character in the framework of a Brownian-motion model for unitary 
matrices. We then study the energy correlations of the conductance and the Wigner time and 
the cross-correlation between these two quantities, making contact with existing calculations 
in the metallic regime. In Sec. ID we undertake similar studies for the localized regime. The 
non-universal parametric correlations of eigenphases and Wigner times can be accounted for, 
at the semi-quantitative level, by a very simple model of resonant transmission through a 
single localized state. We present our conclusions in Sec. [TV|. In Appendix [A], we discuss the 
validity of the assumption on which is based the Brownian motion model, when it is used 
to describe the energy dependence of S. In this case, this amounts to study the underlying 
Wigner-Smith time delay matrix. In Appendix [B] we consider the simple case of one-channel 
scattering, where some exact calculations can be carried out. 



B. Scattering matrix of a disordered region 

We consider an infinite stripe composed of two semi-infinite, perfectly conducting regions 
of width L y (which we define as the leads) connected by a disordered region of same width and 
of longitudinal length L x (see the inset in Fig.|l[a). Assuming non-interacting electrons and 
hard-wall boundary conditions for the transverse part of the wave function, the scattering 
states in the leads at the Fermi energy satisfy the condition k 2 = {nix / 'L y ) 2 +k 2 l , where k is the 
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Fermi wave vector, nn /L y the quantized transverse wave vector, and k n the longitudinal wave 
vector. For a given k, each real transverse momenta labeled by the index n (n = 1, . . . , N) 
define a propagating channel in the leads. Since each channel can carry two waves traveling 
in opposite directions, in regions asymptotically far from the scattering region the wave 
function can be specified by a two 2iV-component vectors, one for each lead (labeled I 
and II). For both vectors, the first (last) N components are the amplitudes of the waves 
propagating to the right (left). In mathematical terms, this reads 

N 1 r 

*i(x,y) = E 7T72 [ A ^ knX + B n e- ik - x ] 4> n {y) (1.1a) 



and 



n=l kr. 



N i 

*n{x,y) = £ [C n e ifc »<*- L «> + D n e-^ x - L ^] cj> n {y). (1.1b) 

n=l kn 



The transverse wave functions are 4> n (y) = y2/L y sin (irny/Ly). The normalization is 
chosen in order to have a unit incoming flux on each channel. The scattering matrix S 
relates the incoming flux to the outgoing flux, 





(1.2) 

With this convention, 5* is a 2N x 2N matrix of the form 

S=(lt)- (1-3) 



The reflection (transmission) matrix r (t) is an N x N matrix whose elements rj, a (tb a ) denote 
the reflected (transmitted) amplitude in channel b when there is a unit flux incident from the 
left in channel a. The amplitudes r' and t' have similar meanings, except that the incident 
flux comes from the right. 

The transmission and reflection amplitudes from a channel a on the left to a channel b 
on the right and left, respectively, are given by [^4 



^Mb) 1/2 / dy' / dy <Pi(y') <f> a {y) G k (L x , y'; 0, y) (1.4a) 



r ba = 5 ab -ih{v a v b )^ 2 / dy' / dy <f>* b (y') <f> a (y) G k (0,y';0,y), (1.4b) 



and 



where v a (v b ) is the longitudinal velocity for the incoming (outgoing) channel a (b). For 
hard- wall boundary conditions, v a = hk a /m, a = a,b. We denote by m the effective 
mass of the electrons. For the transmission (reflection) amplitudes Gk{r f ; r) is the retarded 
Green's function between points r = (x, y) on the left lead and r' = (x', y') on the right 
(left) lead evaluated at the Fermi energy E = h 2 k 2 /2m. Note that, with the convention we 
have taken above, for a perfect, non-disordered sample at zero magnetic field, S is not the 
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identity matrix, but is rather written in terms of transmission submatrices which contain 
pure phases: t ba = t' b * a = S ab exp (ik b L x ). Since we can pass from one convention to the other 



by a fixed unitary transformation, both forms present the same statistical properties (TT[ . 

In our numerical work we obtain the transmission and reflection amplitudes from the 
Green's function of the disordered stripe by a recursive algorithm p5|j26[1 on a tight-binding 
lattice. We typically use a rectangular lattice of 34 x 136 sites and go from the metallic to 
the localized regime by increasing the on-site disorder W. (For details of the simulation see 
Ref. 0.) 

From the transmission amplitudes one can obtain the two-terminal conductance through 
the Landauer formula [27[ 
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Tr (ftt). (1.5) 



Here we adopt units of e 2 /h for the conductance. (Throughout this work we will treat 
spinless electrons and therefore will not include spin-degeneracy factors.) Notice that the 
Landauer formula requires that the sample is a single, complex elastic scatterer. Thus, we 
are ignoring any inelastic process giving rise to a loss of phase coherence. 



C. Eigenphases, conductance, and Wigner time in the metallic and localized regimes 

For a given sample (i.e., impurity configuration) the diagonalization of the scattering 
matrix leads to 2N phase shifts {6{\ as functions of the Fermi energy E. The typical 
dependence is shown in Fig. |l] for the metallic (a) and localized (b) cases for an energy 
interval where new channels are not open (N = 14 in the whole interval). In the metallic 
case the on-site disorder in Anderson units is W = l, yielding a mean-free-path l = 0.2L x , a 
localization length £ = 3L X , and a conductance (thick solid line) which fluctuates around a 
mean value (g)=4.14. 

The localized regime (W = 4, yielding / = 0.02L X and £ = 0.3L X ) exhibits a markedly 
different behavior, with phase shifts showing step-like jumps where the conductance has 
peaks. The peaks indicate that, for certain energies, the probability of traversing the disor- 
dered region is much higher than the average. They are related to the existence of localized 
eigenstates in the sample and transport through the strongly-disordered region is dominated 
by resonant tunneling. 

To illustrate this last point, in Fig. [I] we also show the Wigner time for both metallic 
and localized regimes (thick dashed lines). This characteristic time scale is defined as the 
trace of the Wigner-Smith matrix 



namely, 



7v = Tr(Q). (1.7) 



The unitarity of S trivially implies the fact that Q is Hermitian, and therefore its eigenvalues 
are real. In the case of time-reversal symmetry S T = S, but Q is not necessarily symmetric 
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since in general S' and dS/dE do not commute. Working in a base that diagonalizes S it 
is easy to see that t w admits a simple expression in terms of the energy derivatives of the 
phase shifts, 

The Wigner time can be interpreted as the typical time interval a scattered particle remains 
in the disordered region []20| , p7Tj| . In the localized regime, the Wigner time exhibits the same 
resonant-like structure of the conductance, although peak heights can be relatively different. 
This behavior can be understood in the light of the resonant-tunneling mechanism, since 
each localized state present in the disordered region can trap the electrons for a long time 
||28|| . From this mechanism one can also understand qualitatively why a Wigner-time peak 
can be relatively large when, at the same energy, a conductance peak is small. This may 
happen, for instance, when the tunneling probability rates for channels in lead I is much 
larger than for channels in lead II. We will get back to this discussion in Sec. |TI|. 

A strong correlation between g and r w is also obtained in the metallic regime. The 
correlation in this case is also intuitive, since transport should probe the available density 
of states around the considered Fermi energy. 



II. CORRELATIONS IN THE METALLIC REGIME 



A. Parametric correlations of eigenphases 



In order to characterize the parametric dependence on energy of the set {#/}, we define 
the eigenphase velocity correlator function 



C (AE) = (- 



' d6i(E + AE) ddi{E) 



dE 



dE 




together with the rescaling 



x = AE ^C e (0) 
c e (x)=C e (AE)/C e (0). 



(2.1) 



(2.2a) 
(2.2b) 



The average (• • ■) can be performed over different eigenstates n, over the energy E, or over 
different realizations of disorder. These definitions are analogous to the well-studied case 
]Tj| in which a Hamiltonian H and its eigenvalues {e u } depend on an external parameter <ft 
(say, a magnetic flux) and the eigenenergy velocity correlator is given by 



C £ {A<P) = 



A 2 



'de u (4> + A0) de v {4>) 




(2.3) 



where A denotes the mean level spacing, or inverse density of states around the v-tYt eigen- 
value. For chaotic systems, the universal form of this correlator was checked numerically 
from exact diagonalizations of suitable Hamiltonians [0,in| and after the rescaling 
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x = A(f) ^JC £ (0) (2.4a) 
c £ (x) = C £ (A0)/C £ (O). (2.4b) 

A complete analytical expression for c £ (x) is not available. However, the small and large-x 
asymptotic limits are known exactly from diagrammatic and non-perturbative calculations 



15] , |T6|1 and match accurately the numerical results. In particular, one finds that 



c £ (x) — ► -—. — for x -> oo, (2.5) 

where (3 = 1(2) for spinless systems with preserved (broken) time- reversal symmetry. 

In Fig. H] we present the eigenphase velocity correlation [Eq. (|2.1|) 1 resulting from our 
numerical simulations and the rescaling ( |2.2j ). For weakly-disordered, metallic samples where 
the statistics of the eigenphases at a fixed energy is well described by the Dyson circular 
ensembles [ II[ , we obtain a good agreement with the universal parametric correlation found 



in Hamiltonian systems [pl||T(J. Applying a magnetic field perpendicularly to the stripe, one 
breaks the time-reversal symmetry, causing the parametric correlation behavior to go from 
GOE-like (squares) to GUE-like (circles) |29|. Increasing the disorder (but still remaining 
in the metallic regime) reduces the range of agreement with the universal curve. Further 
increase of the disorder drives the system into the localized regime and away from the 
universal behavior, as we will discuss in the Section III. 

The system-independent form of parametric correlations for energy eigenvalues of random 



Hamiltonians has been studied with a non-linear a model 30 . This treatment has been 



recently extended to show that universality is a property of all systems whose underlying 
classical dynamics is chaotic. For a disordered sample, one finds that universality holds 
when g ^> 1 and the regime is metallic. The same approach has also been used to study the 
statistical fluctuations of the S matrix and the conductance [|1(| of chaotic systems 



under the influence of a generic external parameter. However, there has been no attempt 
to prove analytically the results shown in Fig. [| namely, that the energy correlator of 
eigenphase velocities falls into the analogous curve obtained from energy eigenvalues when 
the number of channels is very large. 

Parallel to the field-theoretical approach, the universality of parametric correlators of 
eigenvalues has also been justified from the hypothesis of a Brownian motion of the eigenen- 



ergies [[33|], with the external parameter playing the role of a fictitious time. The Brownian- 



motion model (BMM) for Hermitian and unitary matrices was introduced by Dyson f53,T3 



and used by Pandey and collaborators |35 in the case of circular ensembles in order to de- 



termine the statistics of the eigenphases at the crossover between different symmetry classes. 
More recently, the BMM has been applied to scattering and transmission matrices describing 
coherent transport through chaotic and disordered systems |36|j37|] . It has been recognized 
in these works that the BMM approach for scattering matrices only obtains correct results 
for a restricted (energy or magnetic field) range, or for sufficiently large number of channels. 
We will illustrate this point in the next subsection, where we consider specifically the energy 
evolution of the phase shifts. 
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B. Brownian- motion model of S and energy-dependent parametric correlators 



The aim of this subsection is to develop a simple description of the parametric statistics 
of phase shifts in the metallic regime numerically investigated (iV^>l). For this purpose, we 
apply a Brownian-motion model to the energy evolution of scattering matrices. As stated 
at the end of the previous subsection, we do not expect to obtain a complete quantitative 
agreement between the BBM predictions and the numerical data - this certainly would 
require a more sophisticated treatment ||. Instead, we only give a justification for the 
agreement observed between the asymptotics of the energy-dependent eigenphase correlation 
functions and analogous curves characteristic of Hamiltonian eigenvalues. 

A unitary matrix S can always be decomposed as the product 

S = Y'Y (2.6) 

of two unitary matrices. In the absence of time-reversal symmetry (/3 = 2), S is just unitary 
and Y and Y' are independent. For time-reversal symmetric systems ((3 = 1), S = S T and 
we have Y' = Y T , where the matrix Y is not unique, but specified up to an orthogonal 
transformation. Any permissible small change in S is then given by 

5S = Y' (iSH) Y, (2.7) 

where 5H is a Hermitian matrix (real symmetric if (3 = 1). This relationship allows us to 
define an invariant measure in the manifold of unitary matrices from the real independent 
components of 5H | 12| , |13| . The isotropic Brownian motion of S occurs when S changes to 



S + SS as some parameter (say, a fictitious time) varies from t to t + St. To construct the 
model, the components of 5H are assumed to be independent random variables behaving 
according to 

= (2.8a) 

and 

f(6H,6H u ) =g/j 5, u , (2.8b) 

with = 1 + 5ij and [i = 1, . . . , 2N + N(2N — 1)(3. As it will be discussed in Appendix [Aj, 
when the evolution of S results from a Fermi energy variation 6E, the Hermitian matrix 5H 
can be given in terms of the Wigner-Smith time delay matrix Q. Thus, the validity of this 
Brownian motion model relies on certain assumptions concerning the statistical properties 
of the various interaction times. A random matrix approach for Q, assuming a maximum 
entropy distribution for a given mean density of eigenvalues, is proposed and some of its 
consequences are numerically checked. Eqs. ( |2.8p are then based on certain simplifications 
which we critically discuss in Appendix [A], at the light of the statistical properties of the 
Wigner-Smith matrix Q. 

The effect of the Brownian motion on the eigenphases of S may be found by second-order 
perturbation theory, 
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56 n 



cot 



Equations fl2.8|) and fl2.9|) then imply that the eigenphases follow the relations 



(2.9) 



(89 n ) = j F{9 n ) St 



and 



(2.10a) 



(59 n 58 m ) = 5 nm — St, 



(2.10b) 



with F(0 n ) = -8W/d9 n and W = - £ n<m In |2 sm[(9 n - 9 m )/2)\. Notice that the coeffi- 
cients / and (3 were introduced to suggest friction and inverse temperature, respectively. 
The eigenphases behave like a classical gas of massless particles on the unit circle, execut- 
ing a Brownian motion under the influence of a Coulomb force F. The joint probability 
distribution of the eigenphases P({9 n };t) follows a Fokker-Planck equation |34| . [13| 



dP 



2iV 



dt f 5 89, 



d fldP 



n=l 



0d9 n 



P 



(2.11) 



For t — > oo we reach an equilibrium situation where the joint distribution becomes that of 
the circular ensemble [TJ, P eq {{9 n }) = C 2N exv[-(3W({9 n })]. 

We now turn to the solution of the Fokker-Planck equation for the joint probability 
distribution of the eigenphases. An exact way to solve Eq. fl2.11|) is to map it into a quantum- 
mechanical Hamiltonian problem of particles in a ring interacting by a two-body potential 
proportional to l/r~, where r 2 c is the length of the chord connecting the pair. This model 



was proposed and studied extensively by Sutherland [38] and yields an exact solution for 
P({9 n };t) [Q. Alternatively, one may use a hydrodynamical approximation, originally 
proposed by Dyson ||40|| , which leads to a non-linear diffusion equation for the average 
density of eigenphases. Here we will adopt this second approach, since it is simpler than 
the first one and sufficient for our purposes. Thus, introducing the t-dependent density of 
eigenphases, 



r 2n 2V 

p{0; t)= d9 x --- d9 2N P({9 l }; t) £ 5(9 - 9 



(2.12) 



and starting from Eq. (|2.11 ), one can derive that 



./ 



dp(e-,t) 

dt 



d_ 
dt 



p(9;t)^J^d9'p(9';t) In 



2 sin 



'9-9 r 



(2.13) 



(This equation is approximate in the sense that it does not take into account accurately short- 
wavelength oscillations in the density.) Since fluctuations in the density are smaller than 
the density itself by a factor of the order 0(1/N), one can linearize this diffusive equation 
by considering small deviations around the homogeneous solution, p(9;t) = p eq + 5p(9;t), 
where p eq = N/tt. One then obtains 
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d5p(9; t) 
di 



d_ 
89 



2n 



d6'D( 



9 



n d 6p(9';t) 
09' ' 



(2.14) 



where the kernel is given by D{9) = —p eq f 1 ln|2sin[(^ — 9')/2]\. The periodicity in the 
eigenphase space calls for a solution of the diffusion equation in Fourier series. Writing 
8p{9- 1) = Et=-oo 5p k (t)e ike , one finds that 



5p k [t) = 5p(0) exp 



7T 



\k\peqt\ 



f J 



(2.15) 



The above result can be readily used in the evaluation of parametric correlators, which will 
be the subject of the rest of this subsection. We therefore introduce the two-point density 
correlator 



R{9,9'-t,t') = (5p(9;t)5p(9';t')) 



(2.16) 



where (• ■ -) eq means that the average is weighted by the joint distribution at the equilibrium, 
Peq{{9i})- Since the average density is constant in both time and angle, there is translation 
invariance in t and 9 at the equilibrium. Consequently, R depends only on the differences 
9 — 9' and t — t'). From Eq. (|2.15 ), we find that the dependence on time and angle can be 
decoupled in A;— space, 



R k (t) = T k exp 



-■n\k\p eq t s 
f , 



(2-17) 



where T k is the Fourier coefficient of the "static" correlator Tiff) = R(9;t = 0), which, for 
circular ensembles, is known exactly ||13|| . In the limit of N 3> 1, one has that T k — \k\/(27i 2 f3) 
for \k\ <ti N, whereas for \k\ 3> iV it saturates at the value T k — N/{2ir) 2 . If we keep 
irpeqt/ f ^> 1/N, we can neglect the contribution of terms with \k\ > N and then easily sum 
the Fourier expansion of R(9; t) to find that 



R(9;t) 



n 2 (3 



Re 



1 



(2.18) 



with z = exp(i9 — 7rp eq t/f). Hereafter we will denote the energy (or time) in terms of a 
dimensionless parameter, with the natural choice being X = E/E s , where E s = ^ f j[j[p eq \ 
The knowledge of R{9\ t) permits the evaluation of other two-point parametric functions. 
Two of them are of particular interest. The first one is the Wigner time correlator 



C T (AE) 



(E + AE)r w (E)) - (r w (E)) 2 

h y ™ \/d9 l (E + AE)d9 m (E)' 
2N 




'd9 l (E)\ ld9 m {E) 



dE 



dE 



dE 



dE 



(2.19) 



which has also been the subject of recent calculations [20,21]. The second parametric func- 
tion is the modified level velocity correlator |16| , |33|| 



2iY 



C(9-9';E-E') = 



'd9 n {E) d9 m {E') 



n,m=l 



dE 



dE' 



5(9 - 9 n (E))5(ff - 9 m (E')) 



(2.20) 



eq 
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It is straightforward to check that 



C T {E) 



71 



d 2 



(2N) 2 dE 2 Jo 



2t 



do 



n) 2 R(9;E) 



(2.21) 



and 



d 2 R d 2 C 



dE 2 de 2 



In Eq. (|2.21| ) one uses that R(8; E) is an even periodic function in 9 and Rk{E) 
k = 0.] Moreover, the two correlators are connected by the relation 



C T (E) 



2tt 



(2Nf 



2- 



al9 C{9;E), 



(2.22) 



for 



(2.23) 



which also implies that C(0; E) = 2{N/n) 2 C T {E). 

Two remarks are pertinent here. First, C T (E) and C(0; E) have asymptotic forms similar 
to Ce(E) in the limit of large E. This is because distinct eigenphases become uncorrelated 
very quickly for large enough energy differences, causing the main contribution to both 
correlators to come from the diagonal terms. The second remark is that, independently 
of the hydrodynamic approximation adopted above, interlevel correlations are completely 
ignored in the BBM [see Eq. ( |2.1(Jb[ )]. This leads to a C T (E) identical, up to a proportionality 
factor, to the one-level velocity correlator Cg(E) for all values of E. Such a coincidence 
between C T (E) and Cg(E) is never observed in practice because interlevel correlations are 
indeed important, particularly among neighboring levels and at small values of E. 

From Eqs. ( p7T7|) , fl2~2l"D , and ( gjgj) one gets 



C T (E) 



n -x 2 



1 + 2X 2 ) 



2(3N 2 E 2 sinh'(X 2 /2) 



(2.24) 



and 



C(9;E) 



C T {E) 



7i 2 {3E\ 



-Re 



z(l - 2X< 



We find the following asymptotic limits: (i) For I < 1, C T (E) 
C(0;E) ~ -A/(tt 2 PE 2 



(2.25) 



-2/((3N 2 E 2 ) and 



for X > 1, C T (E) ~ -AX 2 e- x /((3N 2 E 2 ) and C(0;E) 



—8X e /(tt (3E 2 ). An exponential decay at large distances does not occur for the analo- 



gous correlators involving eigenvalue velocities [33], which actually retains the form (i) for 
arbitrarily large values of the external parameter. As pointed out before |37|j , technically, 
the difference appears because a Fourier series is used to treat phase shifts oscillations (since 
they are bounded by the finite interval [0; 2tt}) instead of Fourier transforms, as in the case 
of energy eigenvalues. Nevertheless, if N^>1, one naturally expects that there should exist a 
range of values of E in which C(9; E) has a form similar to its counterpart for eigenenergies. 
In order to understand this point, let us find the relation between E s and another, more 
commonly used scale, (^(O) -1 / 2 



HUI (see Sec. [flAp . From Eq. ( gTOg) we have that 



11 



C fl (0) 



-1/2 




(2.26) 



which means that E s = Cg(0) 1 ^ 2 y2N/(ir 2 f3). Therefore, if N is large enough, it is possible 
to have E 3> E s and E <C Ce(0)~ 1//2 simultaneously. Consequently, there is an intermediate 
range of values of X where C(0;X) falls exactly into the same power-law, universal curve 
predicted for the case of energy eigenvalues |TE| , |55|] . Deviations in the form of an exponential 
decay occur only for values of X ^ 1, which may be large in units of C*e(0) -1 / 2 / E s . Notice 
that this is still compatible with the restriction E/E s ^> 1/N used to derive Eq. (|2.18f) . 

We stress again that Eqs. (|2.24 ) and (|2.25| ) are approximations. A simple inspection is 
sufficient to prove this point: Because the correlators diverge at E = as 1/E 2 , they cannot 
satisfy the sum rules J °° dE C{6; E) = (for any 9) and / °° dE C T (E) = 0. This limitation 
cannot be overcome within the hydrodynamical approximation we have considered here. In 
fact, Eqs. ( |2.24| ) and ( |2.25[ ) break down for X < (^(O) -1 / 2 . To go beyond this limitation, 
one has to treat Eq. (|2.11|) in a non-perturbative way. 

There is an additional, appealing relation connecting E s to another scale inherent to 



the scattering region. Recall that E s = (h/N)/y /3((r^) — (%) 2 ). The variance of the 
Wigner time was evaluated in Refs. pO|j2T|| using a microscopic formulation to relate S to 
the Hamiltonian of the scattering region, which was modeled as a member of a Gaussian 
ensemble ||. (We remind the reader that Gaussian ensembles are supposed to model the 
statistical properties of ballistic chaotic cavities as well as disordered electronic systems in 
the diffusive regime.) It was shown that 



(V) = M 



\2 



1 + 



P(2N)'< 



(2.27) 



when the scattering region is maximally connected to the external propagating channel (i.e. 
open leads) and N 1. Based on this result, we infer 



E x 



NA 



71 



(2.28) 



This relation says that in the metallic regime of a quasi-lD wire, only one fundamental energy 
scale, the mean level spacing, is required to characterize the decay of energy-parametric 
correlators. Other system-dependent scales, like the Thouless energy E c = TlvfI/L 2 are 
irrelevant. One should recall that, for energy levels, the analogous quantity to C^O) -1 / 2 
is the root-mean square velocity C^O)" 1 / 2 [see Eq. (£§)], which is a direct measure of the 
average dimensionless conductance (g) ~ E c /A when the system is in the metallic regime 
T6| , |45[| . There is no direct relation between C^O) -1 / 2 and C E (0) -1 / 2 and one cannot recover 



any quantitative information about the conductance of the system by calculating Cg(E) 
alone in the metallic regime. 

In Fig. ^ we present the Wigner time correlator (|2.19| ) with the energy rescaled according 
to Eq. (|2.2j ). Notice that the shape of C T is very similar for different symmetry classes ((3 = 1 
and 2). This property also emerges in the context of the stochastic approach to scattering 
when we take the large- iV asymptotics of C T [ 2^, 21 ] : The form of the correlator becomes 
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^ w ' fi(o) [i\"(S ' (2 - 29) 

with T = NA/n. This curve matches reasonably well our data when we let T be a free 
fitting parameter. 

From a semiclassical point-of-view, the existence of time-reversal symmetry in a chaotic 
system does not affect the shape of any energy- dependent, two-point correlator of elements 
of S [41,42| or r w | 43fl . This is because, after energy averaging, these correlators are solely 



determined by the exponential decay with time of the classical probability to escape from 
the scattering region, in which case T _1 is interpreted as the escape rate. The presence or 
not of time-reversal symmetry affects only numerical prefactors in the correlators. The fact 
that large N corresponds to a semiclassical limit becomes clear when we notice that taking 
H — ► for a fixed lead geometry effectively increases the number of propagating channels 

El. 



C. Correlation between the Wigner time and the conductance 

In the previous subsection we discussed the Brownian motion model for S and the impor- 
tance of the probability distribution of the time-delay matrix (further developed in Appendix 
|A|). We now focus on the statistical properties of the trace of Q, the Wigner time t w in the 
metallic regime. As evident from Fig.[l], there are strong correlations between the Wigner 
time and the conductance, that we discuss below. We start by writing the scattering matrix 
in its polar decomposition [|J 

S = UTU. (2.30) 
The 2N x 2N unitary matrices U and U are built out of unitary NxN blocks, namely, 

For systems with broken time-reversal symmetry, the matrices u"> are independent of each 
other. If time-reversal symmetry is preserved, then S is symmetric and one has = u^ T 
and = u^ T . The 2Nx2N matrix T has the block structure 



(2.32) 



T= -n t\ 

\ T Kji 

where 1Z and T are real diagonal NxN whose non-zero elements can be expressed as 

/ A \ 1/2 / 1 \ 1/2 

in terms of the radial parameters A a . The convenience of this representation is that it allows 
the two-probe conductance of Eq. (|1.5| ) to be expressed simply as 
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1 



5 TTV < 2 - 34 > 



0=1 



that is, independent of the unitary matrices U and U. 

In the absence of a magnetic field there is time-reversal symmetry (S = S T , U = U T ), 
and we have 



5 t *s = ut T 

dE 



dE dE dE 



(2.35) 



Since V 2 = I, taking the trace of the above equation simplifies it, yielding 



% ( ,dU dU* T \ 
Tw{E) = -i-^--—V-). (2.36) 

Moreover, because U is unitary, its infinitesimal variations are given by dU = 5U U, where 
SU is antihermitian. Therefore, 

■ & m fSU 5U*\ 
MB ) = TV (--—). (2.37) 

Writing for the block components of U 

du {l) = 5u {l) u {l) and 5u {l) = da {l) + i ds {l) , (2.38) 
I = 1,2, where da® (ds w ) are real antisymmetric (symmetric) NxN matrices, we have 

^' = -sB^ (2 - 39) 



Notice that here (t w ) = 0, since the polar decomposition (|2.30| ) yields the identity for S 



in the absence of disorder. However, with the convention of Eq. ( p..l|) which we took for 
defining our scattering matrix in the numerical simulations, we do not have S = I in the 
absence of disorder. As a consequence, the mean average value of the Wigner time obtained 
numerically is given by the density of states of the disordered region |TT|j46[] . However, as 
expressed in Sec. [I], when we go from one convention to another by multiplying S by a fixed 
unitary matrix, we do not change the statistical properties of the eigenphases. In the same 
way, the constant shift in t w given by the density of states does not change its statistical 
properties. 

One interesting feature of Eqs. (|2.37 ) and ( |2.39| ) is that r w only depends on the infinites- 



imal variations of the unitary matrix U and not on the radial parameters A of the polar 
decomposition. This may seem somehow surprising, given the obvious correlations between 
the Wigner time and the conductance [as shown in Figs, ^.a and [J and the fact that the 
latter depends exclusively on the radial parameters through Eq. (|2.34j ). However, when the 
electron Fermi energy is varied for a given impurity configuration, the corresponding changes 
in A and 8u^ are necessarily correlated since S is constrained by symmetry to move along 
a particular direction in the manifold of allowed polar parameters. In other words, since 
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the Wigner time is given by the energy derivative of S, it must comply with the symme- 
try requirements in the same way as the infinitesimal variations that define the invariant 
measure. 

In Fig. |3] we present the g — r correlator 



(5g(E + AE) 5t w (E)) 



C gr {AE) - //f ^ 2U/2//t _ , 2U/2 , (2.40) 



together with the g—g and t—t correlators. (Here we define St w = (t^) — (t w ) 2 . Analogously 
for Sg.) The correlator of Wigner times has been recently calculated by several authors 
2D| , |2"T[ using supersymmetry methods, while the g — g correlator is only known analytically 



in the metallic perturbative regime 47] when E 3> A. The calculation of the g — r correlator 



will be very interesting, given the numerical evidence provided. 



III. CORRELATIONS IN THE LOCALIZED REGIME 

A. Resonant transmission model for localized transport 

As evident from Fig. p], transport through a localized stripe presents important differences 
for the eigenphases, conductance, and Wigner time as compared to the metallic case of 
Section ||. The peaks of g(E) and r w {E) and the jumps of Qi(E) show that now we are in a 
resonant regime, where the transmission occurs through tunneling into localized eigenstates 
in the bulk of the disordered region. The dependence of transport properties on resonant 
states can be established within the R-matrix formalism ||48| , which allows the scattering 
matrix to be expressed as E9j (see Appendix O) 



W* W 

S nm (E) ~ 5 nm - 2m Y - vv ™ VVm » . (3.1) 

The sum is over the (localized) eigenstates of the disordered region, the matrix elements 
W mu describe the coupling of these states with the different channels in the leads, and 
= 27iJ2n=i \ W nu \ 2 is the resonance total width for the eigenstate v. Equation (|3.1|) is 
valid only to lowest order in F v /A, namely, when resonances do not overlap. (Higher order 
corrections imply traversing the disordered region by sequential tunneling through more than 
one localized state.) Within this approximation, the energy dependence of the conductance 
and the Wigner time appears in the form of Breit- Wigner functions 

(g -iy + n/4 (3 - 2) 

and 

Tw(E) ~ — E {E _ E ^2 + T 2 /r ( 3 - 3 ) 

respectively. The left (r^) and right (r^) partial widths are given by the overlap of the 
corresponding eigenfunction ip u (x,y) with the channel wave functions 4> n (y), i-e., 
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N 



dy 4> n {y)ij u (x = 0,y) 



(3.4) 



and similarly for r^, exchanging x = by x = L. The total width is T u = V® + r£\ The 
coefficients c n = h 2 k n / (2mA) are smooth functions of energy on the scale of A (k n is the 
longitudinal wave vector defined in Sec. |I[B). In the strongly localized case assumed here, 
the typical total width Y u is much smaller than A and only the eigenstate v whose energy 
is the closest to E contributes significantly to the sum. Hence, from now on we will neglect 
any smooth energy dependence in T u and omit the index v. The fitting of g(E) by a single 
Breit-Wigner resonance works quite well for our numerical data when W = 4 or larger; other 
numerical models of disorder [|5(J also yield Breit-Wigner shapes. 

The statistical properties of T are connected to the fluctuations in the eigenfunction 
intensity. Using some very simple arguments one can estimate the probability distribution 
P(r). First we recall that the envelope of a localized state decays as one moves away from 
its center r = (xo,yo)- The scale of this decay is the localization length £, such that 
\4>(r) \ ~ exp (— \r — To|/0- Consequently, we may write that 

r «) ~ e - 2x "^ fff{ij} (3.5a) 

and 

r (r) _ e -2(L-*o)/c (3 5b) 

The factors arise from the fluctuations of the eigenstate on the scale of kp 1 . However, 
since we work with a large number of channels N, Jn follows approximately a Gaussian 
distribution and therefore can be substituted by the its average value. Factorizing away the 
energy scale given by the level spacing A, we define 

r(0 

7(0 = ^A ~ (3 ' 6a) 

and 

7 M = ^1 „ e-^-xo)/^ (3 6b) 

where c is a numerical constant proportional to N. Analogously, we define the dimensionless 
total width 7 = 7O +7 (r) . 

Our statistical assumptions will be the following: 

1. x is uniformly distributed along the disordered stripe, 



P(x ) = \\ (3.7) 



2. xq and ^ are independent random variables; 
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3. z = 2L/£ has a normal distribution [[51 



P(z) = F exp 



2a 2 



(3.8) 



where the mean and the variance are related by a 2 = 2zq. The two first assumptions are 
trivial. The third assumption originates from the standard log-normal distribution of the 

> 5 = < 



dimensionless conductance 51 



with 



var(lng) = —2 (In. 



(3.9) 



The mean value zq is a measure of the disorder, and in the strongly localized regime we 
have Zq^>1. The normalization factor F = [y/nzo(l + ^(a/^o/2))] -1 takes into account the 
fact that z is always positive. (<& denotes the error function.) Nevertheless, to leading order 
in 1/zq, we can ignore this restriction over z and recover the standard Gaussian prefactor 
F ~ [2«/7rzo] -1 . The probability distribution of 7 can be constructed as 



P(i) ~ 2F 



1/2 



ds / dzexp 
Jo 



(z - z f 

4:Za 



d{-f-2e- z/2 cosh(s^)), 



(3.10) 



with s = Xq/L — 1/2. Carrying out one integration, we find that 



, , f z 2 dz 

P(7) ~ 2F / — exp 

J z\ Z 



Vi 2 - 4e" 



(3.11) 



with z\ — 2 In (2/7), z 2 = +00 if < 7 < 1, or z 2 = — In (7 — 1) if 1 < 7 < 2. We cannot 
simplify Eq. (|3.11|) further, and in Fig. ^| we show the result of a numerical integration. 
Obviously, our estimate of P("y) is accurate only when 7<C 1. Working the various asymptotic 
limits, we see that for very small 7 (zx^>Zq), the probability distribution vanishes as 



p( a) oc ^77 exp 



z l 



4:Zn 



(3.12) 



The distribution has a maximum around 7 mp ~ 2e z ^ 2 , which becomes more pronounced for 
increasing disorder, namely, 



P(lm P ) oc — exp 

zq V ^ 

The large values of 7 are not exponentially damped by disorder, as we have 



P(l 



1 

1 oc — . 

z 



(3.13) 



(3.14) 



The knowledge of P(7) allows us to estimate various averages. For instance, let us show 
that the model is consistent. From the Breit-Wigner form of the conductance ( |3.2|) we can 
write 
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fag) 



1 

A 
In 



A/2 
-A/2 



dE In 



7 (0 7 M 



T (0 T W 



(E/cA) 2 + 7 2 /4 



7 2 c 2 



+2 



7c arctan 



7c 



(3.15) 



with the averages taken with respect to -P(7). The first term on the right-hand side gives 
the dominant contribution, 



In 



7 (9 7 (r) 



) = (~z) ^ -z . (3.16) 

Since P{z) selects values of z~zq^>1, the remaining terms of ( 3.15 ) give the next leading- 
order contribution, which is independent of zo, 

2/1 M\ 

+ 2 



In 



7 2 c 2 



+ 1 



21n(2c) + 2. 



The fluctuations of In g are characterized by the correlation function 
CJAE) = (\ng(E + AE) kig(E)) - (\ng(E)) 2 , 



(3.17) 



(3.18) 



with the average running over energy E and disorder. Following the statistical assumptions 
introduced above, we have 



<ln 2 <?> 



1 

A 



A/2 
-A/2 



dE In 2 



7 (0 7 M 



In 2 



7 (0 7 W 



(E/cA) 2 + 7 2 /4 

,2\ \ 



+ tf T. 



+ 2 In 



7 w 7 w 



In ( ^ 



In ( 7 2 c 2 ) 



+ (ln 2 ( 7 2 c 2 ))+4(ln ( 7 2 c 2 

~ z 2 - 2z Q [2 ln(2c) + 1] . (3.19) 

The variance of the distribution of Ing is then given by C g (AE = 0) ~ 2z , showing that our 
resonant model for localized transport is consistent with the standard log-normal distribution 
of the conductance [^TJ characterized by Eq. ( |3.9|) . In subsection |TJ.C, we will apply the 
resonant model to the correlation function (|3.18j ) and will determine its energy-correlation 
length. 



B. Wigner time in the localized regime and correlations with the conductance 



The resonant model can be applied to the Wigner time, whose average is given by 



( T w) 



2N cA 2 



A/2 
-A/2 



dE 



7 



(E/cA) 2 + 7 2 /4 



h 4 / /i N 

( arctan — 

2N A \ V 7 c, 

h 2tt 

2N ~K' 



(3.20) 
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This result also checks the consistency of our model. Notice that although t w fluctuates 
strongly with energy and from sample to sample, (t w ) essentially does not depend on disor- 
der. This is due to the fact that the mean slope of the eigenphases as a function of energy 
is proportional to the density of states [|ll~if46f (or equal to zero, depending on the particular 
convention adopted for S) . Large values of r w at the resonances are compensated by the flat 
parts in between. These large fluctuations are not appropriately represented by the usual 
linear correlation function ( |2.19j ), which, within our resonant model, is given by 

h V 4tt r 2 P(y) 



C ^ = T&h^^T [AE/ (ct A)F + l " M ■ (3 - 21) 

In order to see why, we calculate the second moment of t w : 

, 9 . / % \ 2 87T f°° , P(z) pW dy 

K) = 7T777 TT? / dz 1 



2N J cA 2 Jo ' z 7 7min 7 ^/ 7 2 _ 7 2. n 
k V 2tt 2 



- exp [-z . (3.22) 



\2NJ cA 2 z r V4 

The lower and upper limits of 7 for a given z are 7 m i n = 2e~ z ^ 2 and 7 max = 1 + e~ z - In 
obtaining the above result we have used the fact that the relevant values of z are of the 
order of z ^>l. The (exponentially) large fluctuations of t w reflect a very wide distribution. 
Hence, we will describe these fluctuations in terms of the logarithm of t w . Moreover, for the 
remaining of this section we will adopt the dimensionless Wigner time 

Tw 2N — , (3.23) 



(t w ) t R 

where tu = 2nh/A is the Heisenberg time. Similarly to the calculation of the average 
log- conductance, we write 



1 / r A / 2 
(lnr) =4- ( / dE In 2 



-A/2 



7/(2ttc) 



(E/cA) 2 + 7 2 /4 



~ (In 7) +2 + ln(^\ (3.24) 

The first term in the right-hand side gives the dominant contribution and, to leading order 
in zq, we have 

<lnr)^-|. (3.25) 

Therefore, in the asymptotic limit of z Q ^>l, the mean logarithmic of the conductance and 
the Wigner time are simply proportional to each other, 

(lnr) = \ {\ng). (3.26) 

To investigate the fluctuations around this average value, we proceed in a way analogous to 
Eq. (TO : 
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(In 2 r) 




dE In 2 



7 /(27rc) 



(E/cA) 2 + 7 2 /4 



v [/ 12 
The variance of In r is then given by 

var(lnr) = C T (AE = 0) 



48' 



(3.27) 



(3.28) 



The existence of correlations between the conductance and the Wigner time, as we pre- 
sented in the Introduction (Fig.[l].b), can be easily understood from Eqs. ( |3.2|) -( |3~3"D . More- 
over, our resonant transmission model for localized transport allows us to quantify such 
correlations. For this purpose, instead of working with the cross-correlator (|2.40|) , we define 
the logarithmic correlator 



C gT (AE) 



(5 In g(E + AE) Slxir(E)) 
((51n^) 2 ) 1 /2((5i nr )2)i/2 



(3.29) 



and calculate 

(In glnr) 




dE In 



In 



-A/2 



7 (0 7 (r) 



7 (0 7 M 



(E/cA) 2 + 7 2 /4 



In 



7/(2tt C ) 



In 7) + 



In ('^) +2 + 2h,(r) 



(E/cA)2 + 7 2/4 

111 



7 (0 7 W 



+2[ln(2c) + l](ln 7 ) 



In ( -) + I ln(2c) + 2 

7T/ 2 



(3.30) 



To leading order in zq we have 



C gr (AE = 0) 



(3.31) 



showing that the cross-correlation decays slowly (not exponentially) with disorder. 

In Table I we present our simulations in the localized regime for different disorder (W) 
and number of modes (N) in a stripe with aspect ratio L x /L y = 4. The results displayed 
show a qualitative agreement with the predictions of our resonant model [Eqs. ( ^.16[ ), (fl9|). 
( p.20| ), ( |3.25| ), ( |3.28| ), and ( |3.31|) 1. For the range of disorder that we are able to simulate, 
it is likely that the next-leading order in the large- zq asymptotic expansions of the above 
equations is required for a quantitative agreement. The c— dependent terms should be 
calculated consistently with the value of zq that best describes each sample. We will not 
attempt here this detailed comparison between the model and our numerical simulations 
because we would need a much better statistics than the one we dispose. The relationship 
( p.9| ) between the mean logarithmic conductance and its variance is only approximately 
verified in our numerical simulations. It has already been noticed []50| that in order to 
obtain the agreement with Eq. (|3.9j) one should simulate long stripes with weak disorder. 
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The fluctuations in In g are larger than those in lnr, despite the fact that in our resonant 
model the former depends linearly on z , while the latter goes quadratically with zq. It 
is the difference in the prefactors that makes, for the values of the disorder that we have 
simulated, the fluctuations of In g larger. One interesting aspect of the last column of Table 
I, the value of the cross correlations, is that the results of simulations seem to show an even 
slower dependence on z than that of Eq. (|3.31|) . 



C. Energy and parametric correlations 



The energy-dependent correlation functions for the log-conductance and log-Wigner time 
can be calculated from the resonant model along the same lines that we followed in the previ- 
ous subsections. In particular, the various energy correlation lengths can be estimated from 
the initial curvature of the corresponding correlation functions. For the energy correlator of 
Eq. ( |3.18D we have 



C' g \AE = 0) 




7 (0 T (r) 



-In 




In ( x 2 + 1 



21n(2c) - 1 



(3.32) 



The term proportional to (I/7) gives the dominant contribution and its calculation is anal- 
ogous to that of Eq. ( |3.22j ) , yielding 



C"(AE = 0) 



2tt 2 



cA 2 



~o 



exp -Zq 



(3.33) 



Assuming that the correlation function C g (AE) has approximately a Lorentzian form, its 
correlation energy will be related to the curvature at the origin through 



AE^ 



2C g (0) 



A 



2c 



7T 



Zq exp 



; 2 o 



(3.34) 



For the logarithmic correlation of the Wigner time we have (to leading order in zq) the same 
curvature than for C g (AE) since lnr and \ng have the same energy dependence [Eqs. ( |3.2|) - 
( |3.3D l. The difference in the correlation lengths AE c g and AE^ comes from the different 
values of the variances C g (0) and C T (0), yielding 



AEZ ~ A- 



1 

"4^ V 3 



3/2 



exp 



(3.35) 



Although the resonant model predicts that C T (0) and C g (0) have a quadratic and a linear 
dependence on z , respectively, we find numerically that var(lng) is always larger than 
var(lnr) within the range of disorder and stripe lengths simulated. Therefore, we expect 
that AE c g > AE^, which is indeed the behavior observed in Fig. [5[ The energy-correlation 
length of C gT is intermediate between AE c g and AE°. In all three cases, the correlation length 
shrinks exponentially with Zq, due to the fact that the conductance and Wigner-time peaks 
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become narrower with increasing disorder. This is the reason why we have used in Fig. [5] the 
energy rescaling appropriate for the parametric correlations of eigenphase velocities which 
also turns out to depend exponentially on z (see below). 

The correlation functions in the localized regime can be calculated with the aid of the 
resonant transmission model. The analysis is more involved for eigenphase velocities than for 
g or r w because we do not have a simple expression like (|3.2j ) and ( |3.3|) relating eigenphases to 
energy. However, there are obvious relations between correlators of eigenphases and Wigner 
times. The eigenphase velocity introduced in Eq. ( |2.1|) is simply the diagonal part (I = m) 
of the sum defining the Wigner time correlator of Eq. ( |2.19j ). In particular, Cg(AE) and 
C T (AE) are identical in one trivial limit of the localized case, namely, the N — 1, when we 
know the relation exactly (see Appendix ||): Close to a resonance with energy E u , we have 



0±(E) « ±- + arctan[2(£ 



) 



from which we recover the Breit- Wigner form of 

h r„ 
2 (E-E v y + ri/A- 



>{E) 



(3.36) 



(3.37) 



In the large- N limit and for strongly disordered samples we can assume that the {#/} move 
almost rigidly, repelling each other simultaneously when a resonance occurs. Therefore, there 
is approximately no distinction between diagonal and off-diagonal correlations because all 
eigenphases follow a similar pattern of energy evolution This is confirmed in our numerical 
simulations by the close agreement between C T and Cq over a wide energy range. The large 
fluctuations for the Wigner time translate into very large values of Cg(0), as evident from 
Fig. ^.b. (Notice that the sharp steps of the phase shifts give rise to very large derivatives.) 
The broad distribution of eigenphases forces us to work with the logarithmic correlator 



C e {AE) = (\n 



d9i(E + AE) 
dE 



In 



dME) 
dE 



- ( In 



dOijE) 
dE 



1 \ 2 



but still adopting the rescaling used in the metallic regime [Eq. ( |2.2| )j, 



ce(x 



AE ^C e (0) 
C e (AE)/C e (0). 



(3.38) 



(3.39a) 
(3.39b) 



The average in Eq. ( 3.38 ) is over energy E, disorder, and channel index I. In Fig. |5| we 
present cg(x) for in the localized regime (disorder W = 4). The universality found for the 
corresponding eigenphase velocity of the metallic regime is lost in a way consistent with 
our resonant model. On the other hand, when properly rescaled, we obtain an agreement 
between Ce(x) and the Wigner time correlator. 



IV. CONCLUSIONS 

In this work we have studied parametric correlation functions in disordered quasi-one 
dimensional systems in the metallic and localized regimes. For this purpose we have con- 
sidered the fluctuations in energy of the eigenphases of the scattering matrix, as well as 
fluctuations of the conductance and the Wigner time, and their cross correlation. 
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In the metallic regime, when disorder is weak and the fluctuations of the eigenphases {61} 
are well described by Dyson's circular ensembles [[ll[], the parametric correlations obtained 
from our numerical simulations follow closely the universal behavior discovered by Szafer, 
Altshuler, and Simons |15jl6] for the spectra of chaotic and disordered systems. This finding 
is justified from a Brownian-motion model similar to that developed by Beenakker |33| for 
the energy spectrum. The Brownian-motion model also allows us to obtain the energy 
correlation function of the Wigner time in the large-energy asymptotic limit. We have 
compared these analytical results with our numerical simulations. Our simulations in the 
metallic regime display a strong correlation between the conductance and the Wigner time, 
which arises from the symmetry restrictions to the energy drift of the scattering matrix. 

Transport in the localized regime is resonant-like and, therefore, its statistical properties 
are given by those of localized wave functions in the disordered stripe. From the well-known 
fact that wave functions are (1) localized around centers uniformly distributed in the sample 
(2) with inverse localization lengths following a Gaussian distribution, we recover the basic 
features observed in the simulations where disorder was strong. Our model is consistent 
with the standard log-normal distribution for the conductance and allowed us to estimate 
the typical energy correlation length of conductance fluctuations. The very large fluctuations 
of the Wigner time led us to study the distribution of its logarithm; the variance of this 
distribution was found to be related to that of the conductance. We also investigated the 
correlation between the conductance and the Wigner time as a function of disorder. The 
energy-dependent parametric correlations in the localized regime follow closely the behavior 
of the Wigner time. 

An experimental check of some of the ideas developed in this work has become possible 
with the fabrication of very short insulating wires J52,|53|. In particular, the Breit- Wigner 



line shape characteristic of resonant tunneling has been established and the statistics of 
peak position (as a function of Fermi energy or gate voltage) has been shown to display level 
repulsion ||53|| . This peculiar finding can be been understood within a resonant-tunneling 
picture if the observed peaks correspond to states well connected to the leads, which are 
confined to a small region in center of the stripe and which are therefore not necessarily 
separated from each other by distances much larger than the localization length. But the 
experimental study of short wires poses also the question of whether the single-particle 
approach that we have pursued in this work is valid. At least in a strictly one-dimensional 
sample it is known |54| that two-body potentials drive the system to a non-Fermi liquid 
behavior, leading to a more sensitive dependence on weak disorder [|55|| . Moreover, given 
the recent developments stressing the interplay between disorder and interactions in the 
localized regime it would be interesting to extend the present study in order to 

incorporate the effects of electron-electron interactions in the energy-dependent parametric 
correlations. 
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APPENDIX A: STATISTICAL PROPERTIES OF THE WIGNER-SMITH 

MATRIX 

In this appendix we study the statistical distribution of the Wigner-Smith matrix in 
order to investigate more carefully the applicability of a unitary Brownian-motion model to 
the energy dependence of the scattering matrix. 

For a small energy displacement, the evolution of the scattering matrix can be written 

as 

S(E + SE) = S{E) exp (iSK), (Al) 

where SK is required to be an infinitesimal Hermitian matrix to preserve the unitarity of S. 
To first order in SE, we can write 

SK = -i S^^l 5E = 2 -^Q 5E. (A2) 
dE n 

Consequently, considering now the form ( |2.7| ) of the allowed infinitesimal variations in the 
neighborhood of S, we can identify 

Q5E = ^N Y] {5k) Y > (A3) 

showing that the eigenvalues of Q and SH are proportional to each other. Moreover, if S 
has an isotropic distribution at a given E and remains isotropically distributed for each 
energy in the interval SE, one can see that this property should characterize Q as well. The 
isotropy of S and Q are related. We now try to determine which are the requirements on 
the probability distribution of Q necessary for the validity of the BMM [Eqs. (|2.80 1 at least 
in an approximate sense and for a sufficiently small SE. 

We have investigated the probability distribution of the eigenvalues of Q by a numerical 
simulation on weakly-disordered metallic stripes. Under these conditions, the correlations 



of the eigenphases of S are well approximated by those of the CUE (or COE) ensemble [11 
and the eigenphase velocity correlation agrees well with those found for chaotic Hamiltonian 
systems. In Fig. ^ we present the mean density of eigenvalues of Q and its nearest-neighbor 
spacing distribution (inset) for the cases without and with a time-symmetry breaking mag- 
netic field. Notice that, within our convention, all eigenvalues of Q are positive, consistently 
with their interpretation as typical traversal times through the disorder region. Due to this 
symmetry requirement (that rules out the usual Gaussian ensembles), one of the simplest 
possible random-matrix description of Q is provided by the Laguerre ensembles. The ob- 
served eigenvalue density and the nearest-neighbor distribution (after unfolding M) are not 
incompatible with such an ensemble at first sight. Indeed, one can see a clear level repulsion 
and a good agreement with the Wigner's surmise of the appropriate symmetry class. Notice 
that although Q is not symmetric, the nearest-neighbor distribution in the time-reversal- 
symmetric case is given by the j3 = l Wigner surmise. This is because the submanifold of 
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allowed Q matrices in the presence of time-reversal symmetry can be mapped into that of 
the real symmetric ones by a transformation that leaves the eigenvalues unchanged [B0|. It 



is important to remember that the other set of characteristic times, the energy derivatives 
of the phase shifts {d$i/dE} associated with the eigenchannels of S, do not exhibit nearest- 
neighbor repulsion. Both sets are, however, obviously related since their sum is simply 

T w . 

For a disordered wire, S is not isotropic. (However, the circular ensembles give a good 
description of the phase shift fluctuations in the metallic regime.) Since the isotropy of 5" and 
Q are related, Q cannot be isotropic for a disordered wire. But one can hope that this does 
not matter for the spectral fluctuations of Q, in analogy with what happens for the spectral 
fluctuations of S. This lead us to consider a distribution for the Wigner-Smith matrix that is 
invariant under unitary (orthogonal) transformations, and to propose a simplified maximum- 
entropy ansatz ]5^Q] ; in other words, we adopt a maximum-entropy distribution for Q, 
given the observed mean density of eigenvalues. This will yield a logarithmic interaction 
between eigenvalues and consequently the observed level repulsion. In the usual Coulomb 
gas analogy, one would have a certain (non-parabolic) confining potential for this ansatz 
giving the observed mean density of positive eigenvalues. The problem we face is then quite 
analogous to that of the probability distribution of the radial parameters of the transfer 
matrix |]||. 

It is well known that a maximum-entropy approach with an arbitrary confining potential 
leads to correlations between matrix elements f6lf . The eigenvector isotropy assumption 
allows us to average over the unitary group |62],[13|], yielding 



(Qij) 



(TrQ) 
N 



(A4a) 



and 



(QijQki) 



N 2 



\mn n / \ n mn 



(A4b) 



Above, brackets on the left-hand sides imply averages over the probability distribution of 
Q, while brackets on the right-hand sides indicate averages over the eigenvalue distribution. 
Equation ( |A4a| ) can be trivially set to zero once we suppress the constant drift of the 
eigenphases with energy (or use the convention (t w ) = 0, like in subsection 0.c). The 
non- vanishing correlations of Eq. (|A4b|) are 



(Ql) 



N + 1 



(N-l)( qi q 2 )+2(q 2 ) 



(A5a) 



(Q uQ j j) 



l 



N + l 



N(qiq 2 ) + (q\) 



(A5b) 



(QijQji) 



N+l 



(ql) - (qm) 



* 7^ 3- 



(A5c) 
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It is only in the case of the quadratic confinement potential that we have the GUE re- 
sult (q\q2) = implying a vanishing correlation between different diagonal elements 
[Eq. (|A5bj) 1 and the usual correlation ((QijQji) = (qf)/N) between symmetric matrix el- 
ements [Eqs. (|A5a|) and (|A5cj) 1. Other confining potentials, like the one compatible with 
the observed mean density, result in different two-point correlation functions (giq^) and 
nonvanishing matrix element correlations in all Eqs. (|A5|) . 

If we assume that the eigenvectors of 5H are also isotropically distributed, Eqs. ( |A3| ) and 
( |^5| ) show that there exist correlations between different matrix elements of 5H. Therefore, 
Eqs. ( ^.8| ) cannot be verified exactly and their validity is actually approximate. On the other 
hand, the correlations of the type of Eq. ( |A5b| ) decrease with N and in the large N limit 
the maximum-entropy ensembles provide a local approximation to the Gaussian ones [63] , 
Therefore, it is only in the large N limit that the BMM for S can be appropriate. 

Assuming a nearly Gaussian behavior and comparing Eqs. ( |A5a|) and (|A5cj ) with Eq. 
allows the identifications 

St = {SEf (A6) 

and 

J = f? <TrQ 2 >, (A7) 

of the fictitious time and the friction coefficient of the BMM with the energy and Wigner- 
Smith matrix associated with the scattering process. 



APPENDIX B: SCATTERING IN THE ONE-CHANNEL CASE 

In this appendix we calculate the energy dependence of the eigenphases around a reso- 
nance for the one channel case. The starting point is the S-matrix expression KS| 

S ab (E) = e 2i ^S ab - 2m e ^ + ^ ]T W*^{E)[D~ 1 {E)] tw Whu{E), (Bl) 

with 

D^E) = E5, v - + iir^2W*^(E)W cu (E). (B2) 

c 

W a/1 (E) represents the overlap between the external wave functions (plane waves) and the 
internal eigenfunctions of H (for a time-reversal symmetric system we can choose W aVL to be 
real). The usual approximation is to assume that the energy dependence of the phases f a (E) 
and matrix elements W afl (E) is smooth over the interval where there are many resonances 
(poles) in D~ l (E). Moreover, if |Wa/i| 2 is typically much smaller than the average distance 
between poles, we can expand [D _1 (£')] and perform a unitary transformation in the wave 
functions to obtain 

W* W 

S ab (E) = 5 ab - 27rc£ /9 + 0(r/A), (B3) 
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where T v = 2nJ2c \ W CU \ 2 . Now, specializing for a one- dimensional system (therefore N = 1) 
and looking at energies close to a certain resonance, we can get explicit expressions for the 
coefficients 

1 - i - « ( T- £)/r, < Bte > 
r ' K 1 - i - S%,v, < B4b > 

' K - l-2.(/- W (B4C) 

with a = AkWI v Wrv/T v . Notice that ( [B3| ) keeps S unitary only to lowest order in T/A. 
It is useful to introduce the following parameterization for S: 



r = VR e ir > +x 
r* = VR e v ~ x 
t = iVT e ir > +iK 

t* = iVf e ir >- iK , (B5) 



with R+T = 1. The various quantities appearing above can be determined through Eq. (|B4[) ; 
in particular, 

r]{E) « arctan[2(£ - E v )/T v ] (B6) 

and 

T < E > - i + 4(J-w/ (B7) 

For time- reversal symmetric systems (k — 0), it is easy to find an expression for the eigen- 
phases of S in terms of the new parameters, namely, 

I fr\ 

6± = r] ± arctan \ \ — , mod(7r). (B8) 
Using Eqs. (|B6|) and (|B7|), we then have 

1/2 



9±(E) « arctan[2(£ - £?„)/r„] ± arctan 





a 


2 


1 - 


a 


2 + A(E-E u )yr u 



mod(7r). (B9) 



Notice that, in general, T# ^ Tl and \a\ < 1. This means that the second term on the 
r.h.s of Eq. (|B9| ) varies slower than the first and the difference between eigenphases around 
a resonance is approximately n [Eq. ( |3.36| )]. 
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Note added: While finishing this manuscript, we learnt from C. Beenakker that an exact 
random matrix description of the distribution of the Wigner-Smith matrix Q can be obtained 
in the case that either S remains distributed following the Dyson circular ensembles as E 
varies, or the underlying Hamiltonian is a member of the Gaussian ensembles. This can 
be applied to ballistic chaotic cavities, but not to the disordered wires which we study. 
However, the result is very similar to the maximum-entropy description which we proposed 
in the Appendix [A], notably as far the level repulsion is concerned. The extension of the 
derivation made by Beenakker and co-workers from the ballistic cavity to the disordered 
wire should allow us to see what corrections to the maximum-entropy ansatz proposed in 
Appendix |A] are required. 
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TABLES 



WT 
VV 


AT 




{{ding) ) 




/In -r\ 


lis. ]„ _\2\ 


n 

9 T 


4 


14 


-4.1 


7.0 


1.3 


-0.01 


0.40 


0.67 


4 


10 


-6.2 


9.2 


1.2 


-0.21 


0.62 


0.74 


5 


14 


-9.8 


8.6 


1.2 


-0.43 


0.71 


0.62 


5 


10 


-9.0 


13.8 


1.1 


-0.49 


0.79 


0.60 


6 


14 


-16.1 


29.1 


1.2 


-0.69 


0.80 


0.53 


6 


10 


-17.2 


28.8 


1.3 


-0.99 


0.93 


0.55 


7 


14 


-24.5 


39.3 


1.1 


-0.92 


0.81 


0.52 



TABLE I. Simulations of stripes in the localized regime with aspect ratio L x /L y = 4, mean 
disorder W, and number of channels N. The averages are over energy and impurity configuration. 
(In g) and ((5lng) 2 } are the mean logarithmic conductance and its variance, respectively. The 
mean Wigner time (r) has been normalized to its average value Trh/NA, and therefore the fifth 
column checks the validity of Eq. ( 3.2Cp . (lnr) and ((<51nT) 2 ) are the mean logarithmic Wigner 
time and its variance, respectively. The last column is the value of the cross correlator between 
the conductance and Wigner time calculated as in Eq. (3.29). 
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FIG. 1. Energy dependence of the phase shifts for a quasi-lD disordered stripe (inset) in the 
metallic (a) and localized (b) regimes. The conductance and the Wigner time (thick solid and 
dashed lines, respectively, both in arbitrary units) are smooth functions in the metallic regime and 
exhibit a resonant structure in the localized regime. The energy range in (b) is chosen to facilitate 
the visualization of the resonances. Note the strong correlation between g and t w . 
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FIG. 2. Eigenphase velocity correlator for the metallic regime according to the definitions of 
Eqs. (2.1) and (|2.2| ). The weakly-disordered metallic case (W=l in Anderson units) with (circles) 
and without magnetic field (squares) shows a good agreement with the universal curves character- 
istic of the GUE and GOE, respectively. Increasing the disorder (W=2, no magnetic field, pluses) 
reduces the range of agreement with the universal curve. 
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x = E / C e (0) 1 ' 2 

FIG. 3. Autocorrelations of the conductance g (square) and the Wigner time t w (circle) and 
the cross-correlation between g and t w (triangle) for the metallic case without magnetic field. The 
autocorrelator of Wigner times in the presence of a magnetic field is also plotted (pluses). The 
dotted curve is a plot of Eq. ( |2.29| ) . 
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FIG. 4. Probability distribution for the dimensionless total width 7 for increasing disorder 
(parameterized by zq). Dashed line: zq = 4, corresponding to the case shown in Fig. (|l]).b; 
dash-dotted: zq = 6; dotted: zq = 8; solid line: zq = 10. Inset: blow up of the small 7 region 
showing how the most probable value of the distribution moves towards the origin with increasing 
disorder. 
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FIG. 5. Logarithmic autocorrelations of the conductance (solid), Wigner time (dotted), and 
eigenphase velocity (long-dashed) for the strongly localized case (W=& in Anderson units) without 
magnetic field. The cross-correlation between g and t w is shown by the thick short-dashed curve. 
The use of the logarithmic was required in view of the broad distributions of g and t w (see text) . 
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FIG. 6. Distribution of the eigenvalues of the Wigner-Smith matrix in the weakly-disordered 
metallic regime. The thick and thin lines correspond to the presence or absence of a magnetic field, 
respectively. The inset shows the nearest-neighbor spacing histogram for the same eigenvalues 
(similar convention for the line widths). The Wigner surmises for the GUE and GOE (thick and 
thin dashed lines, respectively) are shown for comparison. 
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